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ABSTRACT 

This  study  introduces  exchange  coefficients  for  wind  stress  (C/>),  latent  heat  flux  (C, ).  and  sensible  heal 
flux  (C's)  over  the  global  ocean.  They  are  obtained  from  the  state-of-the-art  Coupled  Ocean-Atmosphere 
Response  Experiment  (COARE)  bulk  algorithm  (version  3.0).  Using  the  exchange  coefficients  from  this 
bulk  scheme.  Cn.  C, .  and  Cs  arc  then  expressed  as  simple  polynomial  functions  of  air-sea  temperature 
difference  (Tu  -  7\) — -where  air  temperature  (Ta)  is  at  10  m,  wind  speed  (Va)  is  at  10  m,  and  relative 
humidity  (RH)  is  at  the  air-sea  interface — to  parameterize  stability.  The  advantage  of  using  polynomial- 
based  exchange  coefficients  is  that  they  do  not  require  any  iterations  for  stability.  In  addition,  they  agree 
with  results  from  the  COARE  algorithm  but  at  st=5  times  lower  computation  cost,  an  advantage  that  is 
particularly  needed  for  ocean  general  circulation  models  (OGCMs)  and  climate  models  running  at  high 
horizontal  resolution  and  short  time  steps.  The  effects  of  any  water  vapor  flux  in  calculating  the  exchange 
coefficients  are  taken  into  account  in  the  polynomial  functions,  a  feature  that  is  especially  important  at  low 
wind  speeds  (e.g.,  Vu  <  5  ms  ')  because  air-sea  mixing  ratio  difference  can  have  a  major  effect  on  the 
stability,  particularly  in  tropical  regions.  Analyses  of  exchange  coefficients  demonstrate  the  fact  that  water 
vapor  can  have  substantial  impact  on  air-sea  exchange  coefficients  at  low  wind  speeds.  An  example 
application  of  the  exchange  coefficients  from  the  polynomial  approach  is  the  recalculation  of  climatological 
mean  wind  stress  magnitude  from  40-yr  European  Centre  for  Medium-Range  Weather  Forecasts 
(ECMWF)  Re-Analysis  (ERA-40)  data  in  the  North  Pacific  Ocean  over  1979-2002.  Using  ECMWF  10-m 
winds  and  the  authors'  methodology  provides  accurate  surface  stresses  while  largely  eliminating  the  oro- 
graphically  induced  Gibb's  waves  found  in  the  original  ERA-40  surface  wind  stresses.  These  can  have  a 
large  amplitude  near  mountainous  regions  and  can  extend  far  into  the  ocean  interior.  This  study  introduces 
exchange  coefficients  of  air-sea  fluxes,  which  are  applicable  to  the  wide  range  of  conditions  occurring  over 
the  global  ocean,  including  the  air-sea  stability  differences  across  the  Gulf  Stream  and  Kuroshio,  regions 
which  have  been  the  subject  of  many  climate  model  studies.  This  versatility  results  because  C0.  C, .  and  Cs 
arc  determined  for  Vl(  values  of  1  to  40  m  s“‘,  (T„  —  Tv).  intervals  of  —8°  to  7°C,  and  RH  values  of  0%  to 
100%.  Exchange  coefficients  presented  here  are  called  the  Naval  Research  Laboratory  (NRL)  Air-Sea 
Exchange  Coefficients  (NASEC)  and  they  are  suitable  for  a  wide  range  of  air-sea  interaction  studies  and 
model  applications. 


1.  Introduction 

Previously,  Kara  et  al.  (2000,  2002)  introduced  poly¬ 
nomial  functions  for  determining  the  exchange  coeffi¬ 
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cients  of  wind  stress  (CD),  latent  heat  flux  (CL),  and 
sensible  heat  flux  (Cv)  used  in  the  bulk  formulas.  In 
these  algorithms,  the  stability  correction  was  approxi¬ 
mated  by  a  term  depending  on  the  sea-air  temperature 
difference  (Ty  -  Ta),  where  Ta  is  at  10  m,  and  a  second- 
order  polynomial  function  of  the  wind  speed  at  10  m 
above  the  sea  surface  (V(l)  for  Cn  or  the  inverse  wind 
speed  for  CL.  The  algorithms  do  not  require  iteration  to 
account  for  stability;  therefore,  they  are  computation¬ 
ally  efficient  for  use  in  high-resolution  coupled  atmo¬ 
sphere-ocean  circulation  models  and  ocean  general  cir¬ 
culation  models  (OGCMs).  They  are  also  useful  when. 
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every  model  lime  step,  latent  and  sensible  heat  are  cal¬ 
culated  via  the  bulk  formulas  using  both  ocean  model 
sea  surface  temperature  (SST)  and  the  high-frequency 
(e.g.,  6-hourly  or  daily)  wind  speeds.  This  approach 
keeps  model  SST  on  track  without  any  explicit  relax¬ 
ation  to  SST  data  or  a  climatology  (e.g.,  YVallcraft  et  al. 
2003). 

The  algorithms  presented  in  the  two  studies  men¬ 
tioned  above  were  obtained  using  tabulated  exchange 
coefficients,  which  were  published  by  Bunker  (1976) 
and  Isemer  el  al.  ( 1989).  The  use  of  these  two  datasets 
resulted  in  a  few  limitations  in  the  algorithms.  For  ex¬ 
ample,  the  datasets  did  not  provide  realistic  exchange 
coefficients  at  high  and  low  wind  speeds,  predomi¬ 
nantly  to  allow  for  the  biases  inherent  in  estimating 
fluxes  from  ship  observations.  Many  later  studies  indi¬ 
cated  lower  values  for  the  exchange  coefficients  (e.g., 
DeCosmo  el  al.  1996;  Yelland  et  al.  1998,  2002),  in 
general.  As  a  result,  the  algorithms  based  on  these  two 
sets  usually  resulted  in  errors  (5%-25%)  in  estimating 
the  wind  stress  and  latent  and  sensible  heal  fluxes  in 
comparison  to  those  obtained  from  the  standard 
Coupled  Ocean-Atmosphere  Response  Experiment 
(COARE)  version-2. 5b  (v2.5b)  algorithm  (Fairall  el  al. 
1996).  The  COARE  v2.5b  algorithm  determines  the 
turbulent  air-sea  fluxes  through  a  recursive  solution  for 
the  atmospheric  dynamic  stability.  However,  this  algo¬ 
rithm  also  has  some  deficiencies  because  it  was  only 
verified  for  winds  up  to  12  m  s  '.It  is  also  noted  that 
when  Kara  et  al.  (2000)  was  published  and  Kara  el  al. 
(2002)  was  under  review,  the  COARE  v2.5b  algorithm 
was  the  accepted  standard  for  air-sea  flux  estimates  in 
the  peer-reviewed  scientific  literature. 

Fairall  et  al.  (2003)  later  developed  a  new  version  of 
the  COARE  algorithm  (v2.6  in  2001,  renamed  as  v3.0 
in  2003)  with  the  aim  of  providing  more  realistic  ex¬ 
change  coefficients  al  higher  wind  speeds  based  on  the 
data  collected  since  the  release  of  COARE  v2.5b.  The 
new  algorithm  also  takes  additional  observational  data 
(Yelland  el  al.  1998;  Hare  et  al.  1999)  into  account, 
which  suggests  that  a  Charnock  constant  value  of  0.011 
(Smith  1988),  used  in  the  earlier  COARE  v2.5b  algo¬ 
rithm.  is  loo  low  for  higher  wind  speeds.  Thus,  a  linear 
increase  in  the  Charnock  parameter  is  introduced  be¬ 
tween  10  and  18  ms  1  and  a  value  of  0.018  is  used  for 
>18  m  s  1  (Fairall  et  al.  2003).  This  might  overestimate 
the  roughness  at  very  high  wind  speeds  where  Taylor 
and  Yelland  (2001)  predict  a  decrease  in  the  Charnock 
parameter. 

With  the  availability  of  the  improved  COARE  algo¬ 
rithm  v3.0  (Fairall  et  al.  2003),  it  is  now  possible  to 
obtain  more  realistic  Cn  and  CL  values.  This  would 
reduce  the  uncertainties  in  the  datasets  (—  I ()%— 20%) 


used  for  deriving  the  exchange  coefficients  presented  in 
Kara  et  al.  (2000,  2002).  In  addition,  the  use  of  the 
COARE  v3.0  has  the  advantage  of  providing  reliable 
transfer  coefficients  at  low  and  high  wind  speeds.  In  this 
paper,  we  also  provide  an  accurate  formulation  of  air- 
sea  fluxes  using  stability-dependent  exchange  coeffi¬ 
cients,  allowing  the  use  of  ocean  model  sea  surface  tem¬ 
perature  (SST)  in  their  calculation  for  a  wide  range  of 
atmospheric  conditions  over  the  global  ocean.  This  ap¬ 
proach  can  substantially  improve  the  accuracy  of  an 
ocean  model  simulated-SSTs  by  providing  a  physically 
natural  tendency  toward  accurate  SSTs  without  relax¬ 
ation  to  any  SST  data  or  climatology. 

The  main  purposes  of  this  study  are  to  present  (i) 
exchange  coefficients  based  on  the  COARE  v3.0  (sec¬ 
tion  2)  and  (ii)  polynomial  parameterizations  of  CD, 
C/  ,  and  Cs  based  on  these  coefficients  (section  3).  In 
particular,  the  stability  correction  in  these  parameter¬ 
izations  is  approximated  by  a  term  depending  on  the 
(7„  —  TJ  and  a  polynomial  function  of  the  V„  at  10  m 
using  tabulated  exchange  coefficients  as  obtained  from 
the  COARE  algorithm  (v3.0).  A  virtual  temperature- 
based  correction  to  (T„  —  Tx)  is  later  applied  to  the 
polynomial  functions,  thus  including  the  effects  of  wa¬ 
ter  vapor  flux  on  the  exchange  coefficients  (section  4). 
Finally,  exchange  coefficients  presented  in  sections  3 
and  4  are  used  for  calculating  climatological  mean  wind 
stress  magnitude  in  the  North  Pacific  Ocean  over  1 979— 
93  (section  5). 

2,  The  COARE  exchange  coefficients 

The  COARE  algorithm  (v3.0)  presented  in  Fairall  el 
al.  (2003)  is  based  on  previously  published  results  and 
2777  one-hour  covariance  flux  measurements.  Fairall  et 
al.  (2003)  added  4439  new  values  from  field  experi¬ 
ments  between  1997  and  1999  to  test  the  algorithm.  The 
new  observational  dataset  now  dominates  the  database, 
especially  in  the  wind  speed  regime  above  10  ms  ', 
where  the  number  of  observations  increased  from  67  to 
about  800.  After  applying  quality  controls,  the  database 
was  used  to  evaluate  the  algorithm.  The  average  (mean 
and  median)  model  results  agreed  with  the  measure¬ 
ments  to  within  about  5%  for  moisture  from  0  to  20 
ms  For  stress,  the  covariance  measurements  were 
about  10%  higher  than  the  model  at  wind  speeds  over 
15  ms-1  while  inertial-dissipation  measurements 
agreed  closely  at  all  wind  speeds.  Thus,  it  should  be 
emphasized  that  the  COARE  algorithm  (v3.0)  is  con¬ 
sidered  to  be  the  state-of-the-art  algorithm  in  this  pa¬ 
per. 

For  the  purpose  of  this  paper,  the  COARE  algorithm 
(v3.0)  was  modified  to  produce  Cn  and  C,  al  various 
(Ta  —  T,)  and  Vn  intervals  over  the  global  ocean,  fn 
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particular,  the  algorithm  was  run  using  various  values 
of  Tv,  relative  humidity  (RH)  just  above  the  sea  surface, 
T<r  and  Va.  There  is  no  cool  skin  or  warm-layer  calcu¬ 
lation;  thus,  a  skin  sea  surface  temperature  value  is  used 
in  the  algorithm.  The  use  of  a  bulk  value  is  an  approxi¬ 
mation,  which  allows  one  to  perform  the  calculation 
without  specifying  radiative  parameters  or  a  time  his¬ 
tory  of  the  fluxes.  The  roughness  length  for  tempera¬ 
ture  is  assumed  equal  to  that  for  humidity  in  the  algo¬ 
rithm;  thus,  C,  is  equal  to  Cs  in  this  paper. 

In  addition  to  (Ta  —  Tv)  and  Vtn  the  effect  of  any 
water  vapor  flux  through  RH  must  be  taken  into  ac¬ 
count  in  calculating  the  exchange  coefficients  at  low 
winds,  a  topic  which  will  be  discussed  in  section  4  in 
detail.  In  this  section  the  exchange  coefficients  are  ob¬ 
tained  for  a  relatively  low  7S  (10°C)  and  for  humidity 
close  to  or  at  saturation  (RH  =  100%)  to  investigate 
the  impact  of  only  (T„  —  TJ  on  Cn  and  C, .  Note  that 
dividing  the  mixing  ratio  (qn)  by  the  saturation  mixing 
ratio  (qsin)  and  then  multiplying  this  number  by  100 
determines  the  RH  (in  %).  A  7V  value  of  I0°C  was 
chosen  for  consistency  with  the  COARE  algorithm 
testing,  and  because  almost  no  sensitivity  to  the  choice 
of  Ts  was  found.  This  is  true  since  the  only  effect  of  the 
temperature  in  the  COARE  algorithm  enters  through 
the  nonlinearity  of  the  saturation  vapor  pressure  for 
water  vapor,  which  has  some  effect  at  30°C  and  upward. 

Exchange  coefficients  for  Cn  and  CL  are  first  gener¬ 
ated  using  an  RH  value  of  100%  for  various  ( Ta  —  1\) 
ranging  from  -8°  to  7°C:  7°.  6°,  5°.  4°.  3°,  2°.  1°.  0.75°. 
0.5°,  0.25°,  0,  -0.098°,  -0.25°,  -0.75°,  -1°,  -2°,  -3°, 
-4°,  -6°.  — 8°C.  These  values  were  chosen  to  represent 
a  wide  variety  of  ( T(i  -  Ts)  intervals  over  the  global 
ocean.  It  is  noted  that  a  (T„  —  Ts)  value  of  **—0.098° 
gives  the  neutral  stability  value.  The  resolution  of  (Ttt  — 
Ts)  intervals  was  increased  near  the  neutral  stability 
value  to  better  examine  behavior  of  exchange  coeffi¬ 
cients  with  respect  to  stability.  In  addition,  the  10-m 
wind  speed  values  of  l^fO  m  s  1  with  increments  of  l 
m  s  1  were  used  to  determine  coefficients  at  both  low 
and  high  wind  speed  conditions  (e.g.,  tropical  cyclones 
and  hurricanes). 

Figure  1  shows  CD  and  C,  values  obtained  from  the 
COARE  algorithm  when  RH  is  set  to  100%.  Both  Cn 
and  Cf  experience  large  changes  with  wind  speed  at  10 
m  above  sea  surface,  and  they  are  also  quite  different 
with  respect  to  (Ta  -  Ts)  values  for  a  given  wind  speed. 
In  general.  CD  is  always  greater  than  2.0  x  10“"'  at  Va 
>  20  m  s  1  for  any  given  ( Tn  —  7’J  value  (Fig.  la).  On 
the  other  hand,  C,  does  not  have  any  significant  vari¬ 
ability  at  high  wind  speeds  (Fig.  1  b).  In  particular,  C,  is 
almost  constant  with  values  ranging  between  M.2  x 
10"3  and  M.3  x  10  3  for  any  given  (T„  -  7V)  values 


above  20  ms"1.  There  are  strong  variations  in  the  ex¬ 
change  coefficients  with  respect  to  (Ta  -  Ts)  at  very  low 
wind  speeds  (V(l  <  3  m  s_l),  and  this  is  especially  evi¬ 
dent  at  V(l  values  of  l  and  2ms"1  (Figs.  lc,d).  The  Va 
=  I  m  s  1  cases  are  intended  to  show  how  variable  CD 
and  C,  can  be  at  very  low  wind  speed  regimes.  Even 
the  CD  and  C,  values  for  Va  =  I  ms"1  are  quite  dif¬ 
ferent  from  those  for  Vn  -  2  ms-1.  We  actually  ob¬ 
tained  Cn  and  C,  values  from  the  COARE  algorithm 
at  each  1  ms"1  interval  up  to  40  ms-1,  but  they  are 
only  shown  at  2  ms  1  intervals  in  the  figure  for  visual 
clarity.  Although  exchange  coefficients  can  exist  for  Vft 
=  0  ms-1,  they  are  ignored  in  the  analyses  because 
air-sea  fluxes  are  quite  small  under  calm  wind  condi¬ 
tions. 

The  impact  of  stability  through  (Ta  -  Ts)  on  both  CD 
and  C/  is  clearly  evident  from  Fig.  1.  Thus,  exchange 
coefficients  reduced  to  neutral  stability,  based  solely  on 
wind  speed  at  10  m  above  the  sea  surface  (e.g..  Large 
and  Pond  1981),  with  no  stability  dependence  (that  is, 
Ttl  —  =  0  and/or  at  saturated  conditions),  will  not  be 

appropriate  when  they  are  used  in  OGCM  simulations. 
This  becomes  a  critical  issue  in  model  simulations,  es¬ 
pecially  when  studying  the  diurnal  cycle  of  upper  ocean 
quantities. 

Regarding  hurricane-force  winds  (e.g.,  Vtt  >  30 
m  s_l),  there  is  almost  no  direct  experimental  evidence 
for  determining  exchange  coefficients.  There  are  very 
few  direct  observations  of  exchange  coefficients  for 
wind  speeds  greater  than  25  ms  '.As  to  C,  and  Cs, 
there  is  the  problem  of  spray  (e.g.,  Andreas  and  Eman¬ 
uel  2001;  Emanuel  2003).  The  COARE  algorithm  used 
here  includes  an  option  to  define  roughness  length  in 
terms  of  the  sea  state.  It  includes  a  formula  based  on 
wave  age  (Oost  et  al.  2002)  or  one  based  on  significant 
height  and  slope  of  the  waves  (Taylor  and  Yelland 
2001).  The  latter  formula  predicts  significantly  lower 
roughness  at  high  wind  speeds  compared  to  the  former 
one  as  explained  by  Taylor  and  Yelland  (2003)  who 
evaluated  performance  of  the  different  formulas  in  pre¬ 
dicting  the  Cf).  Based  on  evidence  from  laboratory 
studies  (e.g.,  Keller  et  al.  1992)  and  comparisons  with 
field  experiments  (e.g.,  Janssen  1997;  Johnson  et  al. 
1998),  the  formula  proposed  by  Taylor  and  Yelland 
(2001 ),  as  used  in  the  COARE  algorithm,  is  applicable 
over  most  conditions.  It  also  succeeds  at  high  winds  and 
very  short  fetches  in  wind-wave  flumes  when  wave  age 
form u his  fail. 

It  should  be  emphasized  that  there  is  very  little  rea¬ 
son  for  the  exchange  coefficients  obtained  from  the  al¬ 
gorithm  to  change  geographically.  The  main  change 
from  one  area  to  another  would  be  due  to  different 
stability  or  different  temperature  levels,  but  any  such 
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5  10  15  20  25  30  35  40 

Wind  speed  (m  s  l) 


T-Ts 

—  7.00°C 
6.00°C 

—  5.00°C 

—  4.00°C 
3.00°C 
2.00°C 
1.00°C 

—  0.7  5°C 
0.50°C 
0.25°C 
0.00°C 
Neutral 

—  -0.25°C 
-0.75°C 

—  -LOO°C 
- 2.00°C 

~3.00°C 
- 4.00°C 

—  -6.00°C 
- 8.00°C 

6  5  10  15  20  25  30  35  40 

Wind  speed  (ms') 


•  2ms 
4ms 

•  6ms 

•  8ms 
10  m  s’ 
12  m  s 
14  ms' 

•  16  ms 
)8ms 
20  m  s 
22  ms’ 

■  24  m  s 

■  26  m  s 
28  m  s 

■  30  ms 

■  32  m  s' 
34  m  s 

■  36  m  s 

■  38  m  s 
40  m  s 


Air-sea  temperature  difference  (  C) 


Air-sea  temperature  difference  (°C) 


FUi.  I.  Variations  of  exchange  coefficients  with  respect  to  (T„  -  TJ  and  Vtt  values  obtained  from  the  COARE  (v3.0)  algorithm:  (a) 
(  n  values  plotted  against  the  observed  10  m  Vu%  (b)  Cf  values  plotted  against  the  observed  10  m  Vtt\  (c)  C ,,  values  plotted  against  the 
observed  ( /'„  /'J.  where  Tu  is  at  10  m;  and  (d)  C,  values  plotted  against  the  observed  (Tt,  -  7J.  where  T„  is  at  10  m.  7’he  small  panels 

inside  (a)  and  (b)  are  intended  to  show  variability  of  Cn  and  C ,  for  only  V(,  <  10  m  s  respectively.  Dashed  lines  in  (c)  and  (d)  show 
C  n  and  (  ,  values  for  V„  I  m  s  '.  It  is  noted  that  the  bulk  formulas  are  Q,  -  CL  L  p(l  Va  (r/„  —  </,)  for  latent  heat  flux.  =  Cs  Cf, 
()tl  Vtt  -  / \)  for  sensible  heat  flux,  and  7  =  ()l(CnV^  for  wind  stress.  Here.  L  is  the  latent  heat  of  vaporization  (2.5  X  I0h  .1  kg-1). 
(>"  t()0  l\J\  ( 7„  1  273.16)|  (in  kg  m  ').  where  is  the  gas  constant  (287.1  J  kg’ 1  K  ').  and  Ptl  is  the  mean  sea  level  pressure. 


changes  should  he  dealt  with  by  the  algorithm.  The  only 
other  possible  changes  are  1)  the  change  in  gravity, 
which  increases  by  up  to  0.5%  going  from  equator  to 
the  Pole;  and  2)  the  nominal  atmospheric  boundary 
layer  depth  (600  m)  used,  which  is  also  realistic  in  most 
areas.  Another  possible  effect  with  geographic  varia¬ 
tion  is  the  sea  stale.  While  the  effect  of  sea  stale  on 
surface  roughness  is  debatable,  it  is  usually  negligible, 
say  for  Vtt  >  l()ms  1  and  in  deep  water,  as  mentioned 
in  Taylor  and  Ye  Hand  (2001). 


3.  Polynomial  approach  for  exchange  coefficients 

The  exchange  coefficients  (CD  and  Ct)  are  depen¬ 
dent  on  both  Va  and  (Tt,  -  TJ  as  discussed  in  section  2 
(see  Fig.  J).  The  effect  of  RH  on  Cn  and  C,  will  be 
added  to  the  polynomial  functions  as  a  correction  term 
because  including  the  full  stability  conditions  (i.e.,  both 
Ta  -  Ts  and  RH)  in  polynomial  equations  is  not 
straightforward  (see  section  4). 

We  first  investigate  whether  or  not  Cn  and  Cf  can 
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o  COARE  3.0 

—  o.75 c 

—  -0.75  to  0.75  ( ‘ 


Air-sea  temperature  difference  (  C) 


Air-sea  temperature  difference  (  C) 


Air-sea  temperature  difference  (  C) 


Air-sea  temperature  difference  (  C) 


Fig.  2.  Examples  of  applying  optimal  polynomial  fits  to  the  exchange  coefficients:  (a)  Cn  for  V„  ~  1  m  s  (b)  C,  for  V(l  =  I  ms  1 . 
(c)  Cp  for  Vtt  ~  6  m  s  and  (d)  C ,  for  Vtt  -  6  m  s  J.  Note  that  the  exchange  coefficients  shown  with  open  circles  are  obtained  from 
the  COARE  algorithm. 


be  expressed  as  simple  polynomial  functions  of  V tl 
and  (7fl  -  l\)  using  Cn  and  C L  values  when  humidity 
is  close  to  or  at  saturation  (i.e.,  RH  =  100%).  For 
this  purpose.  Cn  and  Cr  values  are  plotted  against 
( 7„  -  7\)  for  each  Vff  value.  As  an  example.  Fig.  2 
shows  variation  of  exchange  coefficients  with  respect 
to  (T(l  -  TJ  for  Vtt  =  I  ms  1  and  Va  =  6  m  s 
suggesting  that  we  can  construct  polynomial  functions 
in  ( 7 n  —  7 s)  with  coefficients  that  vary  with  Vfr  We  split 
the  whole  ( T(l  -  1\)  range  (from  -8.00°  to  7.00°C)  into 


three  subranges.  The  reason  for  this  split  was  that  there 
was  no  simple  (low  order)  polynomial  fit  for  Cn  and  C L 
at  low  wind  speeds  when  using  the  whole  (Ttt  —  Ts) 
range. 

Various  ( 7 a  —  7S)  intervals  were  tested  to  obtain  the 
best  available  lowest-order  polynomial  function  for  the 
exchange  coefficients.  Alter  a  few  tests,  we  found  three 
ranges  for  ( Ta  -  7V):  1)  -8.0()°C  <  7a  -  7,  <  -0.75°C. 
2)  —  0.75°C  <  7„  -  7S  <  0.75°C.  and  3)  0.75°C  <  7„  - 
7\  <  7.D0t,C.  Hereinafter,  for  convenience,  these  inter- 
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vals  will  he  referred  to  as  unstable,  neutral,  and  stable 
cases,  respectively.  As  evident  from  Fig.  2,  it  was  pos¬ 
sible  to  represent  the  neutral  case  using  a  linear  equa¬ 
tion,  and  the  other  two  ranges  (stable  and  unstable 
cases)  by  a  quadratic  equation  for  CD  and  C, .  Obvi¬ 
ously.  a  higher  polynomial  fit  would  fit  belter  for  the 
neutral  case  (e.g..  Figs.  2a, b)  but  the  use  of  the  linear 
function  is  preferred  for  computational  efficiency.  The 
resulting  error  with  respect  to  the  CO  ARE  algorithm 


would  also  be  negligible  because  a  wind  speed  of  1 
m  s'"1  entering  the  bulk  formulations  of  wind  stress  and 
latent  heat  flux  is  also  small  for  these  two  distinct  cases. 

We  obtained  a  single  equation  for  Cn  and  C,  that 
represents  all  wind  speeds  from  Vn  values  of  1-40 
ms-1.  Given  that  a  quadratic  polynomial  is  suitable  for 
unstable  and  stable  cases,  and  a  linear  fit  is  suitable  for 
the  neutral  case,  we  then  write  a  generalized  polyno¬ 
mial  equation  for  C n  as  follows: 


Cn 


'  CIM(  v„)  +  Cm(  V„ )( T„  -  Tt)  +  Cm( Va )( T„  ~  Tf  for 
<  C,AV„)  +  Cm(Vtl)(T„  -  TJ  for 

.  Qj()( K)  +  Crn(V„)(T0  -  n  +  CI)2(VJ(T„  -  Tf  for 


— 8.00°C  <  T„  -  Ts<  —  0.75°C. 
-0.75°C  <  T„  -  T„  <  0.75°C, 
0.75°C  <  T„  -  1\  <  7.00°C. 


(1) 


The  same  type  of  equations  are  also  written  for  C,  as  follows: 


C,  = 


' C,jv„ )  +  Cn(V„)(T„  -  r,)  +  Ci2(VJT„  -  Tf 
'  CIM(V„)  +  CIA( V„)(T„-TS) 

SfJK)  +  CU[V„H T„  -  TJ  +  Cl2{ K„)( 7',,  -  Tf 


for 

for 

for 


-8.00°C  <  Ta  ~  Ts<  -0.75°C, 
-0.75°C  <  T„  -  Ts  <  0.75°C, 
0.75°C  <  T„  -Ts<  7.00°C. 


(2) 


Unlike  Kara  el  al.  (2000,  2002),  there  are  no  limits  set 
on  V (i  in  these  formulations  because  we  were  able  to 
lake  low  and  high  wind  speed  values  into  account  in  the 
algorithms  by  using  V(l  of  I,  2,  and  3ms  1  as  obtained 
from  the  CO  A  RE  algorithm. 

The  next  step  is  to  obtain  C/>(),  C/,,,  and  Cl)2  values 
(and  similarly  C/(),  C, j,  and  C/2  values)  to  be  used  in 
the  generalized  equations  of  ( I )  and  (2).  To  accomplish 
Ih is,  the  example  of  the  curve-fitting  process  shown  in 
Fig.  2  for  Vn  values  of  1  and  6ms"1  is  repeated  for  each 
Va  value  ranging  from  1  to  40  m  s  1  (with  l  ms"1  in¬ 
crements).  Thus,  for  example,  the  quadratic  fits  to  Cn 
for  the  unstable  case  yielded  a  tola)  of  40  Cnih  40  C/>h 
and  40  Cn2  values.  Applying  a  linear  fit  would  be  suf¬ 
ficient  for  unstable  and  stable  cases  when  V(t  >  10 
ms  1 ,  but  for  consistency  a  quadratic  fit  was  used  for 
all  V(l  values. 

The  final  step  is  to  express  C/X),  COJ,  and  CD2  values 
(similarly  for  C/  (),  C,  ,,  and  C/2  values)  as  polynomial 
functions  of  Vn  or  V(l  *,  representing  all  V it  values  from 
I  to  40  ms  However,  there  was  no  single  low-order 
polynomial  equation  for  the  whole  Vtt  range  using  all  40 
of  the  C„0,  C/; , ,  and  C/)2  values.  This  is  not  suprising 
because  both  C n  and  Ct  vary  significantly  at  low  wind 
speeds  as  shown  before  (see  Figs.  Ic,d).  Therefore,  we 
used  two  V(l  intervals:  I)  V(f  <  5  ms"1  and  2)  Va  >  5 
ms  '.  Figure  3  shows  C/X)  and  C/{)  values  for  the  Va  ^ 
5  ms  1  case.  While  a  cubic  fit  worked  best  for  C/)0,  a 
quadratic  lit  is  sufficient  for  CM).  In  the  unstable  case, 
we  can  then  write  Cn(]  =  1.89100  -  0.71820  Va  + 


0.19750  Vl  -  0.01790  V*  and  C/  ()  =  2.07700  -  0.39330 
V„  +  0.03971  Vl 

Polynomial  coefficients  for  CD  and  CL  are  given  in 
Table  1  for  unstable,  neutral,  and  stable  cases.  Note 
that  the  best  fit  is  sometimes  obtained  when  using  V(l  1 
rather  than  V(,  as  the  dependent  variable.  Since  Va  1  is 
being  used  in  some  cases,  there  must  be  a  minimum  for 
Vtr  Here  the  minimum  V(l  considered  is  I  ms"1.  It  is 
also  important  to  note  that  we  obtained  all  polynomial 
functions  such  that  they  match  up  at  Va  =  5  m  s  l,  and 
at  (Ta  -  Ts)  values  of  -0.75°  and  0.75°C.  as  well. 

4.  Relative  humidity  effects  on  exchange 

coefficients 

The  polynomial  functions  presented  in  section  3  were 
determined  using  exchange  coefficients  from  the 
COARE  v3.0  algorithm,  which  was  run  for  various  (Tit 
—  Ts)  values  at  RH  =  100%,  ignoring  the  impact  of  RH 
on  the  exchange  coefficients.  This  simply  implies  that  in 
our  approximations  so  far,  by  far  the  major  effect  of 
water  vapor  on  CD  and  CL  is  its  temperature-depen¬ 
dent  impact  on  stability,  which  is  very  large  in  the  Trop¬ 
ics  and  decreases  with  decreasing  absolute  tempera¬ 
ture.  Because  of  the  exponential  increase  of  saturation 
vapor  pressure  with  temperature,  in  regions  such  as  the 
Tropics,  humidity  has  a  first-order  effect  on  the  stabil¬ 
ity,  and  since  the  Tropics  tend  to  be  light  wind  areas,  on 
the  fluxes  as  well. 

Here,  our  purpose  is  to  include  full  atmospheric  sta- 
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Wind  speed  at  10  m  (m  s  ])  Wind  speed  at  10  m  (m  s  ') 

Fig.  3.  The  wind  speed  dependence  values  of  Cou  and  C,u.  Results  are  shown  for  V()  <  5  ms'1  cases  for  each 

given  ( T„  -  7\)  interval. 


bility  in  the  polynomial  functions  (see  section  3)  by 
using  water  vapor  in  addition  to  (Ta  -  TJ  and  Vtr 
The  exchange  coefficients  calculated  using  these  full 
stability  conditions  will  be  correct  over  a  wide  range 
of  conditions.  For  example,  our  polynomial  approach 
is  applicable  to  conditions  where  the  water  vapor  (in¬ 
cluding  RH)  makes  a  significant  contribution  to  the 
buoyancy  flux,  that  is,  when  the  full  impact  of  stability 
on  the  air-sea  duxes  is  included.  Similarly,  one  would 
not  expect  any  significant  error  when  Tu  is  warmer 
than  Tv,  but  the  water  vapor  results  in  unstable  strati¬ 
fication. 

The  polynomial  functions  for  C D  and  C,  shown  in 
section  3  assume  saturated  conditions  (i.e.,  RH  = 
100%).  However,  water  vapor  effects  can  be  taken  into 
account  via  a  virtual  temperature-based  correction  to 
(Tit  -  Ts).  The  virtual  air  temperature  (Tlxt)  is  calcu¬ 
lated  from  the  mixing  ratio  at  10  m  above  the  sea  sur¬ 
face  (c/J  using  Tw  =  Tit  [1  +  0.61  qa(Ta)\  with  Ta  in 
Kelvins  (K).  One  possible  approach  would  be  recast  the 
polynomial  functions  for  Cn  and  C,  to  be  functions  of 
(TiXt  -  Tus)  rather  than  (T„  —  Tv).  Since  we  already  have 
polynomial  functions  for  saturated  conditions,  we  in¬ 
stead  apply  a  correction  that  is  zero  when  RH  =  100%: 

(T(l  ~  TX  -(T{t~  Ta)  ~  0.61  7>/sal  ~  qj,  (3) 


where  c/SM  is  the  saturation  mixing  ratio  at  Ta  from  a 
simplified  version  of  the  original  formulation  for  satu¬ 
rated  vapor-pressure  presented  by  Lowe  (1977),  which 
has  a  computational  advantage  over  Buck  (1981). 

The  polynomial  functions  for  Cn  and  C,  applied  to 
(T„  —  Ts)t  are  obviously  identical  to  the  original  ap¬ 
proach  or  saturated  conditions  and  are  in  good  agree¬ 
ment  with  the  COARE  algorithm  (v3.0)  for  all  relative 
humidities  (not  shown).  Once  the  correction  is  applied 
to  the  polynomial  functions,  exchange  coefficients  arc 
calculated,  ft  is  emphasized  that  in  the  polynomial  func¬ 
tions  one  could  use  very  fine  increments  for  (Ta  —  7\) 
rather  than  those  used  in  the  original  COARE  algo¬ 
rithm  (see  section  2)  to  calculate  CD  and  C,  because  all 
polynomials  are  continuous  for  any  given  ( Tn  —  TJ,  Vtr 
and  RH. 

Using  the  tabulated  exchange  coefficients  obtained 
from  the  polynomial  functions,  we  now  investigate  the 
impact  of  RH  (0%—  1 00%)  in  calculating  CD  and  Cr  it 
is  clear  that  assuming  an  RH  value  of  100%  results  in 
significant  errors  for  Va  =  1  ms-1  (Fig.  4).  The  impact 
of  RH  on  CD  and  C,  is  clearly  seen  when  the  difference 
between  T(l  and  Ts  is  small  (e.g.,  -2°C  ^  T(,  -  Ts  < 
2°C).  Such  a  result  clearly  demonstrates  that  atmo¬ 
spheric  stability  can  be  predicted  incorrectly  when  ig- 
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Taiii  i  I.  Optimal  polynomial  cocfficienls  for  Cn  and  CL.  They  are  shown  for  V(t  <5ms  1  and  V(l  >  5  m  s  1  separalely.  Note  that 
Cn  and  C,  coefficients  are  continuous  at  Vtf  =  5  ms  (Tit  -  Ts)  =  0.75°C.  and  ( Ttr  -  T\)  =  ().75°C. 


Polynomial  coefficients 


Wind  speed 

Cubic 

0 

1 

2 

3 

-8.00°  C  s  T„  - 

Ts  <  0.75°C 

Va  £  5  m  s  1 

(  /») 

Vtl 

1.891  00 

-0.718  20 

0.197  50 

-0.017  90 

(  />) 

Vu  1 

-0.006  30 

-0.302  80 

0.312  00 

-0.121  00 

C»2 

Vtl  1 

0.000  44 

-0.017  69 

0.013  03 

-0.003  39 

V(l  . '  Sms  1 

(  no 

v„ 

0.649  70 

0.069  93 

3 .54  c -5 

-3.43e-6 

Cm 

vtt  1 

0.003  83 

-0.275  60 

-1.091  00 

4.946  00 

Cn: 

v,;' 

— 4.83e-5 

0.007  71 

-0.255  50 

0.765  40 

-0.75°C  <  T„  - 

■  7\  £  0.75°C 

Vu  <  5  m  s  1 

Cno 

v„ 

0.977  40 

-0.256  60 

0.104  80 

-0.010  97 

Cm 

V(l  1 

0.205  10 

-1.903  00 

1.133  00 

-0.265  80 

Vu  -*  5  m  s  1 

(  no 

Vti 

0.543  80 

0.083  16 

-0.000  49 

3.09e-6 

Cm 

vtt  1 

-0.016  69 

0.573  80 

- 1  2.240  00 

32.530  00 

().75°C  <  Tn  - 

T  <  7.00°C 

Vtl  v  5  m  s  ' 

(  />() 

vu 

-0.066  95 

0.313  30 

-0.001  47 

-0.004  06 

Cm 

v„ 1 

0.099  66 

-2.116  00 

4.626  00 

-2.680  00 

Cn: 

V(l  1 

-0.024  77 

0.272  60 

-0.555  80 

0.313  90 

Vt,  ~~  5  m  s  1 

Cno 

V" 

0.558  10 

0.081  74 

-0.000  45 

2.67c-6 

Cm 

t;1 

-0.005  59 

0.209  60 

-8.634  00 

18.630  00 

Cm 

Vn  1 

0.000  60 

-0.026  29 

0.212  10 

0.775  50 

Quadr. 

-8.00°C  <  Ta- 

T  <  -0.75°C 

V„  £  5  m  s  1 

(  to 

V„ 

2.077  00 

-0.393  30 

0.039  71 

c,  1 

Va 

-0.289  90 

0.073  50 

-0.006  27 

(7.2 

v,, 

-0.019  54 

0.005  48 

-0.000  49 

V„  ^  5  ill  s  1 

(7.U 

vtl 

1.074  00 

0.005  58 

5.26e-5 

(7.i 

Vtl  1 

0.006  91 

-0.224  40 

- 1 .027  00 

0.2 

vo  1 

0.000  19 

-0.002  18 

-0.101  00 

-0.75°C  <  T(l  - 

7\  <  0.75°C 

Vn  £  5  ni  s  1 

(/,, 

vtt 

0.858  00 

0.097  43 

-0.010  56 

(7 1 

Vu 

-1.927  00 

0.734  50 

-0.077  06 

V  t  ^  5  m  s  1 

(  7  n 

vu 

1 .023  00 

0.009  61 
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-0.29250 

0.549  80 
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0.210  30 

-5.329  00 

(7.2 

V"  1 

0.00.1  55 

-0.062  28 

0.509  40 

noring  RH  dependence  (i.e.,  assuming  humidity  close 
to  or  at  saturation)  in  calculating  CD  and  C,  at  low 
wind  speeds.  Overall,  the  direct  effect  of  decreasing  RH 
is  to  generally  increase  exchange  coefficients  as  seen 
from  the  scalterplols  of  Cn  and  CL  for  V a  =  i  ms  1 
(Fig.  5).  For  example,  if  RH  were  20%  and  one  as¬ 
sumed  an  RH  value  of  100%,  then  CD  would  be  under¬ 
estimated  with  =6  times  less  than  its  actual  value  (0.176 
x  10  3  instead  of  0.762  X  10  3)  for  Vtl  =  1ms  1  and 
Ta  “  Ts  =  l.()"C  (Table  2).  Similarly,  the  underestimate 
of  Cf  would  be  ^7  times  less  than  its  actual  value 
(0.220  x  10  ’  instead  of  1.626  X  10  for  Va  =  I  ms"1 
and  7’,  -  I\  =  0.7°C  (Table  2).  On  the  other  hand,  it 


must  be  noted  that  the  effect  of  RH  on  exchange  coef¬ 
ficients  is  relatively  small  for  high  wind  speeds  as  can  be 
seen  in  CD  and  C,  values  for  V (l  -  6  m  s  1  (Table  3). 
While  the  water  vapor  effect  is  most  important  at  low 
wind  speed,  the  reader  is  cautioned  that  low  wind  speed 
implies  small  wind  stress  and  small  latent  and  sensible 
heat  fluxes.  If  the  wind  stress  or  latent  and  sensible  heat 
flux  values  are  small  to  start  with,  even  an  extremely 
large  change  in  their  strength  might  still  mean  small 
values  after  the  RH  effect  is  taken  into  account.  For  ex¬ 
ample,  latent  heat  flux  at  Vtt  =  1  m  s  1  is^l.9  W  m  2 
(-=1 1.8  W  m  2)  when  RH  =  100%  (RH  -  20%).  In  both 
estimates  we  use  Tu  -  Ts  =  2°C,  and  q(,  -  q „  =  3  g  kg- '. 
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Air-sea  temperature  difference  (MC)  Air-sea  temperature  difference  (  C) 


Flu  4.  Exchange  coefficients  plotted  against  various  ( Tn  -  7J  values  at  V„  —  \  ms  1  as  determined  from  the 
polynomial  functions:  (a)  Cr)  and  (h)  C, .  This  figure  demonstrates  variability  of  exchange  coefficients  with  water 
vapor.  Results  are  shown  for  RH  values  of  20%.  40%,  60%.  80%.  and  100%. 


5.  Wind  stress  magnitude  over  the  North  Pacific 
Ocean 

As  an  example,  results  from  calculating  wind  slress 
magnitude  based  on  the  stability  dependent  Cn  is  dem¬ 
onstrated  in  the  North  Pacific  Ocean  (Fig.  6).  All  analy¬ 
ses  are  based  on  6-hourly  meteorological  data  from  an 
archived  product,  the  40-yr  European  Centre  for  Me¬ 
dium-Range  Weather  Forecasts  (ECMWF)  Re-Analy¬ 
sis  (ERA-40;  Simmons  and  Gibson  2000).  A  climato¬ 
logical  mean  wind  stress  magnitude  is  formed  over 
1979-93.  This  time  period  was  chosen  to  be  consistent 


with  ERA- 1 5.  which  has  been  extensively  used  in  pre¬ 
vious  ocean  modeling  studies  (c.g.,  Hurlburt  and  Hogan 
2000;  Kara  et  al.  2004).  In  addition,  all  fields  shown  are 
interpolated  onto  the  1/12°  Hybrid  Coordinate  Ocean 
Model  (HYCOM)  domain  to  demonstrate  possible 
problems  when  using  atmospheric  variables  (e.g.,  wind 
stress  fields)  from  coarse-resolution  archived  products 
(e.g.,  E125°ERA-40  data)  in  forcing  much  finer 
OGCMs  (e.g.,  1/12°  HYCOM).  Basic  features  and  a 
description  of  HYCOM  are  given  in  Kara  et  al.  (2005). 

We  first  form  a  climatological  mean  total  wind  stress 
magnitude  field  (Fig.  6a)  using  zonal  and  meridional 
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given  RH  values  against  those  for  RFI  =  100%. 
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Tabu-  2.  Exchange  coefficients  lor  varying  RH  values  and  (Ttl  —  7\)  intervals  at  V(,  -  1  m  s  1  as  determined  from  the  polynomial 
junctions.  Also  given  are  Cn  (RH)/Cy,  (RH  =  100%)  ratios.  In  the  table,  for  example.  Cn  (80% )/('/>  (RH  100%)  is  the  ratio  of  Cn 
for  RH  80%  to  that  for  RH  =  100%.  Note  that  all  exchange  coefficients  must  be  multiplied  by  10 
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ro 

Cn  for  varying  RH 

ai  V„  = 

1ms  1 
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(RH)/C„ 
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0.032 
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0.039 

1.00 
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LOO 
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1.00 

LOO 

wind  stress  components  that  are  directly  output  at 
6- hourly  lime  intervals  from  ERA-40.  This  field  is  ob¬ 
tained  from  the  planetary  boundary  layer  model  in  the 
ECMWF  atmospheric  model.  There  is  excessive  noise 
in  the  wind  stress  magnitude  near  the  coastal  bound¬ 
aries,  such  as  along  the  eastern  Pacific  and  in  the  Sea  of 
Japan.  This  noise  is  a  major  motivation  for  recalculat¬ 
ing  the  wind  stress  magnitude  rather  than  using  the  one 
output  from  ERA-40.  Thus,  our  purpose  is  to  eliminate 


the  noise  shown  in  Fig.  6a,  so  that  this  field  can  be 
suitable  as  an  atmospheric  forcing  field  for  an  OGCM 
simulation.  We  recompute  the  wind  stresses  |r  —  (rv, 
rv) ]  with  a  bulk  formulation  using  CD  values  as  derived 
in  sections  3  and  4  as  follows: 

%  =  p(tCnu{ir  +  ir)1,2,  (4) 

tv  -  p(tCnvitr  +  ir)1  (5) 

t  -  p(lCn{ir  +  ir).  (6) 
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Tabu-:  3.  Exchange  coefficients  for  varying  RH  values  and  ( Tt(  -  7'J  intervals  at  V tt  —  6  m  s  1  as  determined  from  the  polynomial 
functions.  Also  given  are  Cn  (RH)/C„  (RH  100%)  ratios.  In  the  table,  for  example,  Cn  (80% )/C’„  (RH  =  100%)  is  the  ratio  of  Cn 
for  RH  -  80%  to  that  for  RH  =  100%.  Note  that  all  exchange  coefficients  must  be  multiplied  by  10  ' 
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1.12 

1.16 

1.20 

3.000 

0.781 

0.809 

0.838 

0.868 

0.899 

0.931 

1.00 

1.04 

1.07 

1.11 

1.15 

1.19 

2.000 

0.870 

0.899 

0.929 

0.960 

0.992 

1.007 

1 .00 

1 .03 

1.07 

1.10 

1,14 

1.16 

1 .000 

0.970 

0.997 

1.012 

1 .052 

1.113 

1.140 

1 .00 

1.03 

1.04 

1 .08 

1.15 

1.18 

0.700 

0.997 

1.013 

1.053 

1.1 12 

1.140 

1.152 

1.00 

1.02 

1 .06 

1.12 

1.14 

1.16 

0.500 

1 .006 

1 .039 

1 .095 

1.133 

1.149 

L 1 58 

LOO 

1.03 

1.09 

1.13 

1.14 

1.J5 

0.200 

1 .043 

1.100 

1.134 

1.150 

1.158 

1. 1 73 

1 .00 

1.05 

1 .09 

1.10 

1.1 1 

1.12 

0.000 

1 .084 

1.128 

1.147 

1.154 

LI  68 

1.182 

LOO 

1.04 

1.06 

1 .06 

1.08 

1.09 

-0.098 

1.108 

1.136 

1.150 

1.159 

1.173 

1.187 

1 .00 

1.03 

1.04 

1 .05 

1.06 

1.07 

-0.200 

1.121 

1.143 

1.152 

1.164 

1.178 

1.191 

1 .00 

1.02 

1.03 

1.04 

1 .05 

1.06 

-0.700 

1.152 

1.163 

1.176 

1.J89 

1.201 

1.213 

1 .00 

1.01 

1.02 

1 .03 

1.04 

1.05 

- 1 .000 

1.165 

1.178 

1.191 

1 .203 

1.215 

1.226 

1.00 

1.01 

1.02 

1.03 

1.04 

1 .05 

-2.000 

1.215 

1.226 

1 .236 

1.246 

1 .256 

1 .265 

1.00 

1.01 

1.02 

1 .03 

1.03 

1.04 

-3.000 

1 .260 

1.268 

1.277 

1.285 

1 .292 

1 .300 

1 .00 

1.01 

1.01 

1 .02 

1.02 

1.03 

-4.000 

1.298 

1.305 

1.31  1 

1.318 

1.324 

1.330 

1.00 

1.01 

1.01 

1.02 

1.02 

1.02 

6.000 

1 .356 

1 .360 

1.364 

1.367 

1.371 

1.374 

1.00 

1.00 

1.01 

1.01 

1.01 

1.01 

8.000 

1.39! 

1.391 

1.391 

1.391 

1.391 

1.391 

1 .00 

1.00 

LOO 

LOO 

l  .00 

1 .00 

where  V(l  =  ( ir  +  ir)l/2.  r  -  (r;  +  17.) I/2.  and  it  (v)  is  the 
zonal  (meridional)  wind  speed  component  (ms  '). 
Density  of  the  air  at  the  air-sea  interface  is  obtained 
using  p„  -  100  PJ\ (/„  4-  273.I6)|  (in  kg  m  '). 
where  is  the  gas  constant  (287.1  J  kg  1  K  ')  and 
P(f  is  the  mean  sea  level  pressure.  All  parameters 
needed  for  calculation  of  Cf).  p„,  and  wind  stress  mag¬ 
nitude  (i.e.,  Ttr  Ts.  qtr  //,  u,  Pn)  are  obtained  from 
ERA-40  as  mentioned  earlier.  Wind  stress  magnitude 


is  computed  at  each  6-h  interval  and  averaged  over 
1979-93. 

The  unrealistic  noise  along  the  coastlines  is  greatly 
cleaned  up  by  using  the  bulk  formulation  with  stability- 
dependent  Cn  based  solely  on  Ta  -  Ts  (Fig.  6b).  This  is 
also  true  when  calculating  Cn  depending  on  both  Ttt  - 
1\  and  water  vapor  (Fig.  6c).  The  spatial  patterns  and 
approximate  magnitude  of  the  mean  wind  stress  mag¬ 
nitude  from  the  bulk  formulation  usually  agree  with 
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(a)  Wind  stress  magnitude  from  ERA40  (N  m  2) 


18  588885388  si 

(b)  Wind  stress  magnitude  without  RH  (N  m-2) 


(d)  Wind  speed  at  10  rn  from  wind  stress  (ms  l) 


(e)  Ratio  of  wind  stress  magnitude:  (c)/(b) 
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(e)  Wind  stress  magnitude  with  RH  (N  m  2)  (f)  land-sea  mask  used  in  ERA40 


Fie;.  6.  Climatological  mean  wind  stress  magnitude  over  1979-93  (a)  as  obtained  directly  from  the  ERA-40  output  itself,  (b)  as 
computed  from  live  bulk  formulation  using  Cn  ignoring  RH  effects  (see  section  3),  arid  (c)  as  computed  from  the  bulk  formulation  using 
C  n  with  the  inclusion  of  RH  effects  (see  section  4).  Also  shown  are  (d)  wind  speed  at  10  m  above  the  sea  surface  calculated  from  wind 
stress  field  shown  in  (c).  (e)  the  ratio  of  wind  stress  magnitude  calculated  using  (c)/(b).  and  (f)  the  land-sca  mask  from  ERA-40,  and 
note  that  a  grid  point  is  either  land  (I )  or  sea  (0).  giving  values  between  0  and  I  for  interpolated  the  land- sea  mask.  All  fields  shown 
are  interpolated  on  the  1/12°  Mercator  grid  used  in  the  Hybrid  Coordinate  Ocean  Model  (HYCOM). 
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those  from  ERA-40  itself.  Basin-averaged  mean  wind 
stress  magnitudes  are  0.058,  0.034,  and  0.037  N  m  2  for 
Fig.  6a,  Fig.  6b,  and  Fig.  6c,  respectively.  The  relatively 
large  basin-averaged  mean  value  of  0.058  N  m  2  in  Fig. 
6a  is  almost  twice  as  large  as  those  in  Figs.  6b  and  6c. 
This  is  due  generally  to  the  unrealistically  large  wind 
stress  magnitudes  near  the  land-sea  boundaries,  result¬ 
ing  from  excessive  noise  in  the  original  ERA-40  prod¬ 
uct. 

The  recomputed  wind  stress  magnitude  can  also  be 
used  for  determining  Va.  Thus,  excessive  noise  near  the 
coastal  boundaries  can  be  removed  from  the  Vn  field  as 
well  (Fig,  6d).  Here,  assuming  neutral  Cn  values,  Vit  is 
approximated  using  Vm  =  l/[p„  X  10  3  (0.926  +  3.2 18r 
-  5.23 It2  +  3.236r )]  for  r  <  0.771 1  and  Vu  =  l /[p„  x 
10  3  (1.461  +  0.485t  -  0.092 r2  +  0.007r3)|  for  t  > 
0.7711.  The  use  of  the  wind  stress  magnitude  field  in 
Fig.  6c  instead  of  that  in  Fig.  6b  did  not  make  any 
significant  difference  in  computing  the  V(l  because  the 
ratio  of  these  two  fields  is  close  to  1  over  the  most  of 
North  Pacific  Ocean  (Fig.  6e),  revealing  that  the  impact 
of  water  vapor  in  calculating  Cn  is  small  on  the  clima¬ 
tological  time  scales.  However,  there  are  differences 
when  Vtt  is  small,  and  this  is  consistent  with  the  results 
shown  in  section  4.  The  reader  should  be  cautioned  that 
there  could  be  impacts  of  water  vapor  on  Cn  on  shorter 
time  scales  (e.g.,  daily  or  shorter),  and  this  is  especially 
true  at  low  wind  speeds  and  low  relative  humidities  (see 
section  4). 

Finally,  Fig.  6f  shows  the  land-sea  mask  used  in  the 
original  1.125°  X  1.125°  Gaussian  grid  for  ERA-40  in¬ 
terpolated  onto  the  1/12°  HYCOM  domain.  For  ex¬ 
ample.  a  contour  value  of  0.6  implies  that  the  wind 
stress  magnitudes  on  the  model  grid  are  ^60%  con¬ 
taminated  by  the  wind  stress  magnitude  over  land. 
Some  grids  near  the  coast  are  represented  as  land 
points  in  ERA-40  but  most  of  them  should  be  ocean. 
Such  problems  in  atmospheric  forcing  fields  were  also 
noted  other  studies  (e.g.,  Metzger  and  Hurlburt  2001; 
Kara  et  al.  2005).  Thus,  it  is  clear  that  unrealistic  wind 
stress  magnitudes  near  such  coastal  region  errors  (see 
Fig.  6a)  are  largely  due  to  the  land-sea  mask  and  oro- 
graphically  generated  Gibb’s  waves.  Gibb’s  waves  are 
caused  by  the  spectral  truncation  of  the  orography  to 
triangular  truncation  (P.  Kallberg,  ECMWF,  2003.  per¬ 
sonal  communication).  When  averaging  fields  origi¬ 
nally  created  in  the  ERA-40’s  spectral  geometry,  traces 
of  the  originally  minute  Gibb's  waves  become  more 
visible.  For  example,  Metzger  (2003)  noted  such  atmo¬ 
spheric  forcing  problems  in  an  ocean  model  study  of  the 
South  China  Sea.  On  the  other  hand,  the  atmospheric 
variables  (e.g.,  /'„.  //,  u,  etc)  in  the  ERA-40  product  are 
from  post  processed  fields  where  the  spectral  traces 


have  been  removed,  allowing  one  to  calculate,  for  ex¬ 
ample,  wind  stress  magnitude  from  these  variables  as 
we  did  in  the  North  Pacific  Ocean  (Figs.  6b,c). 

6.  Summary  and  conclusions 

Ocean  models  or  coupled  atmosphere-ocean  models 
can  be  directly  forced  with  net  surface  heat  fluxes  ob¬ 
tained  from  archived  weather  products.  However,  re¬ 
calculating  air-sea  fluxes  that  are  dependent  upon  the 
model  SST  is  typically  sufficient  to  keep  model- 
simulated  SSTs  on  track  without  any  explicit  relaxation 
to  SST  data  or  a  climatology  (although  assimilation  of 
SST  data  can  further  increase  the  model  SST  accuracy). 
This  process  requires  reliable  and  efficient  parameter- 
izations  of  stability-dependent  exchange  coefficients 
over  the  global  ocean,  a  topic  that  is  investigated  in  this 
study. 

We  present  simple  and  computationally  efficient  pa- 
rameterizations  of  exchange  coefficients  for  wind  stress 
(C n)  and  latent  heat  flux  (C/ ),  which  are  of  particular 
value  not  only  for  high-resolution  global  OGCM  appli¬ 
cations  but  also  for  various  ocean  mixed-layer  and  air- 
sea  interaction  studies.  The  polynomial  functions  for 
exchange  coefficients  are  determined  after  running  the 
COARE  v3.0  algorithm  for  wind  speed  magnitudes  at 
10  m  above  the  sea  surface  ranging  from  1  to  40  m  s  l. 
The  stability  dependence  in  the  exchange  coefficients 
was  accomplished  by  using  air-sea  temperature  differ¬ 
ence  and  air-sea  mixing  ratio  difference  through  RH 
for  a  given  wind  speed  al  10  m  above  the  sea  surface.  It 
is  shown  that  ignoring  water  vapor  effects  in  the  calcu¬ 
lation  of  exchange  coefficients  can  result  in  serious  er¬ 
rors  at  low  wind  speeds,  generally  a  serious  underesti¬ 
mation  of  the  exchange  coefficients  for  V„  <  5  m  s  '. 
As  an  application  of  our  methodology,  climatological 
mean  wind  stress  magnitude  is  calculated  using 
6-hourly  meteorological  variables  from  ERA-40  in  the 
North  Pacific  Ocean  over  1979-93.  Such  an  analysis  can 
identify  possible  problems  existing  in  atmospheric  forc¬ 
ing  parameters  (e.g.,  wind  stress  magnitude),  which  are 
directly  output  from  coarse-resolution  archived 
weather  products  (e.g..  ERA-40)  and  demonstrates 
how  these  problems  can  be  reduced  using  our  method¬ 
ology.  The  use  of  exchange  coefficients  from  the  poly¬ 
nomial  functions  in  the  wind  stress  bulk  formulation 
not  only  reproduces  wind  stress  magnitudes  approxi¬ 
mately  comparable  to  the  original  ERA-40  field  in  the 
interior  of  the  North  Pacific  Ocean,  but  also  eliminates 
unrealistic  and  excessive  noise  near  the  land-sea 
boundaries  that  exist  in  the  original  ERA-40  field. 

Finally,  the  motivation  for  this  paper  was  to  present 
sufficiently  accurate  transfer  coefficients  for  air-sea 
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fluxes  l hat  were  also  computationally  efficient  enough 
to  he  used  for  climate  model  applications.  Although  the 
earlier  COARE  algorithm  (v2.5b)  required  ^20  itera¬ 
tions,  the  latest  version  (v3.0)  iterates  3  times  or  less. 
However,  even  having  three  iterations  at  each  time  step 
would  result  in  significant  computational  inefficiency 
for  a  fine  resolution  global  OGCM  that  runs  many 
years  (e.g.,  decadal  to  inlerdecadal  time  scales).  Ex¬ 
change  coefficients  based  on  our  derived  polynomial 
functions  can  be  calculated  at  times  lower  compu¬ 
tational  cost  than  the  more  involved  computations  of 
COARE  algorithm  (v3.0)  with  three  iterations  for  at¬ 
mospheric  stability.  We  also  note  here  that  the  sensible 
heal  flux  coefficient  Cs  is  assumed  to  be  equal  to  C,  in 
this  paper  because,  for  simplicity,  the  roughness  length 
for  temperature  is  assumed  equal  to  that  for  humidity 
in  the  algorithm.  Further  investigation  concerning  vari¬ 
able  roughness  lengths  is  beyond  the  scope  of  this 
study. 

The  exchange  coefficients  for  air-sea  fluxes  pre¬ 
sented  in  ihis  paper  are  called  the  Naval  Research 
Laboratory  (NRL)  Air-Sea  Exchange  Coefficients 
(NAS EC).  Two  simple  FORTRAN  programs  (one  for 
C />  and  one  for  C, )  along  with  the  tabulated  exchange 
coefficients  based  on  various  wind  speed,  relative  hu¬ 
midity  and  air-sea  temperature  difference  intervals  can 
be  downloaded  from  the  Web  site  (available  online  at 
http://wvvw7320.nrlssc.navy.mil/nasec/nasec.html). 
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